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Abstract — We consider the two user wireless Multiple Access 
Relay Channel (MARC), in which nodes A and B want to 
transmit messages to a destination node D with the help of a relay 
node R, For the MARC, Wang and Giannakis proposed a Com- 
plex Field Network Coding (CFNC) scheme. As an alternative, 
we propose a scheme based on Physical layer Network Coding 
(PNC), which has so far been studied widely only in the context 
of two-way relaying. For the proposed PNC scheme, transmission 
takes place in two phases: (i) Phase 1 during which A and B 
simultaneously transmit and, R and D receive, (ii) Phase 2 during 
which A, B and R simultaneously transmit to D, At the end of 
Phase 1, R decodes the messages xa of A and xb of B, and 
during Phase 2 transmits /(xa^xb), where / is many-to-one. 
Communication protocols in which the relay node decodes are 
prone to loss of diversity order, due to error propagation from 
the relay node. To counter this, we propose a novel decoder which 
takes into account the possibility of an error event at R, without 
having any knowledge about the links from A to R and B to R, 
It is shown that if certain parameters are chosen properly and if 
the map / satisfies a condition called exclusive law, the proposed 
decoder offers the maximum diversity order of two. Also, it is 
shown that for a proper choice of the parameters, the proposed 
decoder admits fast decoding, with the same decoding complexity 
order as that of the CFNC scheme. Simulation results indicate 
that the proposed PNC scheme performs better than the CFNC 
scheme. 

I. Background and Preliminaries 

We consider the two user Multiple Access Relay Channel 
(MARC) shown in Fig. [T] Source nodes A and B want to 
transmit messages to the destination node D with the help of 
the relay node R. All the nodes are assumed to have half- 
duplex constraint, i.e., the nodes cannot transmit and receive 
simultaneously in the same frequency band. In addition to the 
presence of direct link, communication paths exist from the 
source nodes A and B to the destination node via the relay 
node R. As a result, in a two user MARC channel, a diversity 
order of two can be achieved, if the transmission scheme is 
chosen properly. 

A. Background 

In a wireless network, due to the superposition nature of 
the wireless channel, signals interfere at the nodes. Avoiding 
this interference by making the nodes transmit in orthogonal 
time/frequency slots incurs a loss of spectral efficiency. The 
concept of physical layer network coding, in which the nodes 
are allowed to transmit simultaneously resulting in interfer- 
ence, was first introduced in 1 1 1. Physical layer Network Cod- 
ing (PNC) has been shown to outperform traditional schemes 
which involve orthogonal transmissions yj- ||4||. So far, most 



of the works on physical layer network coding have mainly 
focussed only on the two-way relay channel. In this paper, we 
propose a scheme based on PNC for the MARC. 

In a two-way relay channel, in order to ensure unique 
decodability at the end nodes, the network coding maps used 
at the relay node should satisfy a condition called the exclusive 
law (5). These maps satisfying the exclusive law form a 
mathematical structure called Latin Squares and the properties 
of Latin Squares have been used to obtain the network coding 
maps in a two-way relay channel |[6J- |[8|. An interesting 
connection between the proposed PNC scheme for the MARC 
and the two-way relay channel is that the network coding map 
used at R needs to satisfy the exclusive law, for the proposed 
PNC scheme for the MARC to achieve a maximum diversity 
order two. 

Wireless relay networks, in which the relay nodes decode 
the messages are prone to loss of diversity order, due to the 
forwarding of erroneous messages. Various methods have been 
proposed in the literature to avoid this loss of diversity order. 
Cyclic Redundancy Check bits are used so that the nodes 
forward only those packets which are decoded correctly |[9|. 
Some works assume the knowledge of all the instantaneous 
fade coefficients or error probabilities associated with the 
intermediate nodes at the destination node, with the decoder 
at the destination using this knowledge to ensure full diversity 
fTO| , plj . Another method used widely is to use a scaling 
factor at the relay nodes which depends on the fade coef- 
ficients, with the scaling factor indicated to the destination 
using pilot symbols 1 12], 1 13 |. The proposed scheme does not 
suffer from the disadvantages of any of the above methods, yet 
ensures maximum diversity order. This is achieved by means 
of an efficient choice of the transmission scheme and a novel 
decoder used at the destination D. 

A Complex Field Network Coding (CFNC) scheme for the 
MARC was proposed in |13|. In this paper, as an alternative, 
we propose a PNC scheme. As observed in p3| , when R 
transmits a many-to-one function of A's and B's transmission 
during the relaying phase and minimum squared Euclidean 
distance decoder employed at D, a. loss of diversity order 
results. In this paper, we show that for the proposed PNC 
scheme, making the source nodes also transmit during the 
relaying phase, combined with a novel decoder which is not 
minimum squared Euclidean distance decoder, ensures the 
maximum possible diversity order of two. Furthermore, if 
certain parameters are chosen properly, the proposed decoder 
for the PNC scheme can be implemented with a decoding 



complexity order same as that of the CFNC scheme. 
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Fig. 1: The two user MuMple Access Relay Channel 

The following are the main advantages of the proposed PNC 
scheme over the CFNC scheme proposed in |T3l: 

• In the CFNC scheme, R transmits a complex linear 
combination of A's and B's messages and the signal 
set used at R during the relaying phase has points, 
where M is the size of the signal set used at A and 
B. In contrast, since the proposed PNC scheme uses a 
many-to-one map, the signal set used during the relaying 
phase has only M points and hence the PNC scheme is 
expected to perform better than the CFNC scheme, which 
is confirmed by the simulation results. 

• In the CFNC scheme, R uses a scaling factor which 
is a function of the fade coefficients of the A-R and 
B-R links, which needs to be indicated to D using 
pilot symbols. Since the proposed PNC scheme does not 
involve any such scaling factor, there is no need of such 
pilot symbols. 

Notations: Throughout, vectors are denoted by bold lower case 
letters and matrices are denoted by bold capital letters. The set 
of complex numbers is denoted by C. CA/'(0,cr^) denotes a 
circularly symmetric complex Gaussian random variable with 
mean zero and variance and A/'(0, a^) denotes a Gaussian 
random variable with variance cr^ . For a matrix A, A'^ and A* 
denotes its transpose and conjugate transpose respectively. For 
a matrix A, rank (A) denotes its rank and det(A) denotes its 
determinant. For a complex number x, x* denotes its conjugate 
and \x\ denotes its absolute value. For a vector v, ||v|| denotes 
its Euclidean norm. The total transmission energy of all the 
three nodes is assumed to be equal to Eg and all the additive 
noises are assumed to have a variance equal to 1. By SNR, 
we denote the transmission energy Eg. For a signal set <S, AS 
denotes the difference signal set of 5, AS = {x X I X ^ X ^ 
S}. The all zero matrix of size n x n is denoted by On- IE(X) 
denotes the expectation of X. 



B. Signal Model 

Throughout, a quasi-static fading scenario is assumed with 
the channel state information available only at the receivers. 

Let S denote the signal set of unit energy used at A and 
B, with M = 2^ points, A being a positive integer. Assume 



that A and B want to transmit A-bit binary tuples to D. Let 
/i : F2 5 denote the mapping from bits to complex symbols 
used at A and B. 

For the proposed PNC scheme, transmission occurs in 
two phases: Phase 1 during which A and B simultaneously 
transmit and, R and D receive, followed by the Phase 2 during 
which A, B and R transmit to D. 

Phase 1: Let xa = I~i{sa), xb = I^{sb) ^ S denote 
the complex symbols A and B want to convey to D, where 
sa^sb ^ 1^2- During Phase 1, A and B transmit scaled 
versions of xa and xb respectively. The received signal at 
R and D during Phase 1 are respectively given by. 



Vr = hAR\/EsaxA + hBR\/EshxB + zr, and 

HDi = hADV^CiXA + JlBD^/E^bXB + , 



(1) 



where a, 6 G C are constants and the additive noises zr 
and are assumed to be CA/'(0, 1). The fade coefficients 
are Rayleigh distributed, with Jiar ^ CA/'(0, a^^), Hbr ^ 
CA/'(0,a|^), Had - CM({),(j\b) and Hbd - CA/'(0,a|^). 

Let G 5^ denote the Maximum Likelihood (ML) 

estimate of [xaiXb) at R based on the received complex 
number yR, i.e.. 



[xa.Xb) = arg 



mm \yR 



hAnVElaxA - hBR\/E~shxB 



Phase 2: During Phase 2, A and B transmit scaled versions 



of Xa and xb respectively and R transmits xr = f{: 



%R /^R 
A^^B 



), 



where / : 5^ ^ 5 is a many-to-one map. The received signal 
at D during Phase 2 is given by, 

= hAuVE^cxA + hBoVE^dxE + hRD\^XR + zd2- (2) 

where c, d G C are constants and the additive noise is 
assumed to be CA/'(0, 1). The fade coefficient Hrd is assumed 
to be CA/'(0,cr|^). 

In order to ensure that the total transmission energy at the 
nodes A and B is equal to Eg^ the constants a, 6, c and d are 
chosen such that |ap + |cp = 1 and |6p + |dp = 1. 

From ([T]) and (|2]), the received complex numbers at D during 
the two phases can be written in vector form as. 



IAD 



f^RDl 



YD 



aXA CXa 
bXB dXB 
Xr 



+ [zd^ ZD2] ' (3) 

V ' 

ZD 

The matrix C{xa,xb^xr) in ^ is referred to as the code- 
word matrix. The restriction of C{xa^xb^xr) to the first 
two rows, denoted by Cr{x AjXb) is referred as the re- 

axA CXa 
bxB dxB 
a Ax A cAxa] , 
bAxs dAxsl ' "^^''^ 
Ax A J Ax B G AS are referred to as the restricted codeword 
difference matrices. 



stricted codeword matrix, i.e. Cr{xA,XB) 
The matrices Cr{AxAj Axb) 



From ([3]), the vector can also be written as, 
Yd = V^^^hWA + V^^BhWB + y/E~sXRh.W^ + zd, 



where the matrices Wa 
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, Wb = 
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and 



W 



R 



are referred to as the weight matrices at node 





lo ij 
A, B and R respectively. 

The contributions and organization of the paper are as 
follows: A novel decoder for the proposed PNC scheme is 
presented in Section II A. In Section II B, it is shown that 
the decoder presented in Section II A achieves a maximum 
diversity of two if and only if the following two conditions 
are satisfied: (i) the map / satisfies the so called exclusive 
law and (ii) the constants a, 6, c and d are such that the 
restricted codeword difference matrices have full rank for all 
non-zero values of Ax a and Ax b • In Section III, the condition 
under which the proposed decoder admits fast decoding is 
obtained. It is shown that when the weight matrices Wa 
and Wr (or Wb and Wr) are Hurwitz-Radon orthogonal, 
the proposed decoder admits fast decoding, with the decoding 
complexity order same as that of the CFNC scheme proposed 
in fT3| . Simulation results which show that the proposed PNC 
scheme performs better than the CFNC scheme are presented 
in Section IV. 

II. A Novel Decoder for the Proposed PNC Scheme 
AND ITS Diversity Analysis 

In Section II A, a novel decoder for the proposed PNC 
scheme is presented. In Section II B, the condition under which 
the proposed decoder offers a maximum diversity order two 
is obtained. 

A. A Novel Decoder for the Proposed PNC Scheme 

Consider the case when D uses the minimum squared 
Euclidean distance decoder given by, 

( ^rj, ) = arg min \\yD^ — hAD\/E~sa xa — hBD\/E~sh xb\^ 

+ \yD2 - hAD\/E~sC XA - hBD\/E~sd Xb - h rd \/E~s f {x A , X b)^^ 

Since the above decoder does not consider the possibility of 
error events at the relay node, it does not offer maximum 
transmit diversity order two. 

Alternatively, we propose a novel decoder which considers 
the possibility of error events at R, given by. 



The relay transmits a wrong network-coded symbol, inde- 
pendent of {xa-,xb)-, if the joint ML estimate at the relay 
is such that xr = f{x^^x^) ^ /{xa^xb)- Under 
this condition, the optimal ML decision metric at D is given 
by m2{xA,XB)' The relay transmits a wrong network-coded 
symbol with a probability which is proportional to at 
high SNR. Hence to the metric m2{xA-,XB)-, a correction 
factor \og{SNR) is added and the minimum of mi{xA-,XB) 
and \og{SNR) + m2{xA,XB) is taken to be the decoding 
metric at D. 

The CFNC scheme proposed in (13] uses minimum squared 
Euclidean distance decoder, which has a decoding complexity 
of 0{M'^). Since the decoder given in ^ involves minimiza- 
tion over three variables xa^xb and xr, it appears as though 
the decoding complexity order is 0{M^). In Section III, it 
is shown that by properly choosing the constants c and 

the decoding complexity order can be reduced to (9(M^) 
which is the same as that of the CFNC scheme. 

The diversity analysis of the decoder given in ^ is pre- 
sented in the next subsection. 

B. Diversity Analysis of the Proposed Decoder 

The following theorem gives the condition under which 
the proposed decoder for the PNC scheme offers maximum 
diversity order two. 

Theorem 1: For the proposed PNC scheme, the decoder 
given in ^ offers maximum diversity order two if and only 
if the following two conditions are satisfied: 

1) The map / satisfies the condition called exclusive law 
given by. 



arg min {mi{xA,XB), 

ixA,XB)^S'^ 



\og{SNR)^m2{xA,XB)}, (4) 



where the metrics mi {xa^xb) and 1712 (xajXb) are given in 
([5]) and ^ respectively, at the top of the next page. 

The idea behind the choice of this decoder is as follows: 
If the relay transmits the correct network-coded symbol, the 
optimal ML decoding metric at D is given by mi{xA-,XB)- 



(8) 



f(xA,XB) / f(xA,XB), for XA / x^, V G <S, 
f(xA,XB) / fixA^x'^), for XB / x'^ , V G <S. 

2) The restricted codeword difference matrices 
Cr{AxA,AxB) have full rank, WAxa / 0, Axb 7^ 
and Axa^Axb ^ AS. 
Proof: See Appendix. ■ 

Note that condition 2) does not demand full rank for all the 
restricted codeword difference matrices. In fact, whatever may 
be the choice of a, 6, c and it is impossible to obtain full rank 
for the restricted codeword difference matrices of the form 
Cr(AxA,0) (and also Cr(0,AxB)), since both the entries 
of the second row of Cr(Ax^,0) are zeros. It suffices to 
ensure full rank for only those restricted code word difference 
matrices for which both Ax a and Axb are non-zeros. 

It is easy to find a map / satisfying the exclusive law, since 
all maps satisfying the exclusive law form Latin Squares |[6|. 

Definition 1: [14] A Latin Square L of order M on the 
symbols from the set = {0, 1, • • • , t — 1} is an M x M 
array, in which each cell contains one symbol and each symbol 
occurs at most once in each row and column. 
Two simple examples of Latin Squares are the Modulo-M 
Latin square and the Bit-wise XOR Latin Square. In the 
Modulo-M Latin Square, a cell in the M x M is filled with 
the modulo-M sum of the row index and column index. In the 
bit- wise XOR Latin Square, a cell in the M x M array is filled 



11 {XA, Xb) = \yDi — hAD \/E~sa XA — hBD \/E~sh Xb I + \yD2 — hAD \f^sC XA — flBD \/E~sd Xb " ^fLD y^E^f^XA, xb)\ , 



1712 (x A, xb) = lyoi — hAoV^sa X A — hBoV^sb xb\ + min s \yD2 — ^ad V ^sC c 

\ ' ^R9^f{xA^^B)'^R^'S ^ 



/ — / — Pi 

XA — hBDV Esd Xb — hRDV EsXr\ > 



ms {xa.xb) = \yDi - hADvEsa xa 



■hBDvEsbxBl + min 
I x^es 



^^yD2—^AD\/E^CXA — hBDV'^E^dXB—hRDV'^E^XR^ |. 



(5) 

(6) 
(7) 



the bit-wise exclusive OR of the row index and column index 
represented in binary, after binary-to-decimal conversion. For 
M = 4, the Modulo-4 Latin Square and the Bit-wise XOR 
Latin Square are as shown in Fig. [2] 
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(a) Modulo-4 Latin 
Square 



(b) Bit-wise XOR 
Latin Square 



Cr(AxA, Axb) is 



Fig. 2: Examples of Latin Squares of order 4 

In the following example, a choice of a, 6, c and d is pro- 
vided, which ensures that the restricted codeword difference 
matrices have full rank for all non-zero values of Ax a and 
Axb. 

Example 1: Choosing a = l,6=^,c = and d = 

the restricted codeword difference matrices are of the form 
AxA 
_-^^AxB j=,AxB_ 
full rank for all Axa^Axb / since det {Cr{AxA^ Axb)) = 
-^AxaAxb ^ 0. 

A sufficient condition under which the decoder given in ^ 
admits fast decoding is obtained in the next section. 

in. A Fast Decoding Algorithm for the Proposed 
Decoder 

In this section, it is shown that by properly choosing the 
constants a, 6, c and the decoder given in ^ can be 
implemented efficiently by means of a fast decoding algorithm. 

Before presenting the algorithm, we introduce some nota- 
tions. 

The points in the signal set S are denoted by 5^ , 1 < z < M. 
From ([3]), the vector y^^ can be written as. 



Yd 



ahAD 
cHad 



bhBD 
dhBD 







Xa 
Xb 
xr 



ZD 



He 



The matrix Hgq can be decomposed using QR decomposition 
as Heq = QR, where Q is a 2 x 2 unitary matrix and R = 
[Ri is a 2 X 3 matrix, with Ri being upper- triangular of 
size 2x2 and Y2 being a column vector of length 2. Let r^j 
denote the {i^jY^ entry of R. 

Define yo = Q^Yd^ = [vdi yD2V - 



Also, let 

(I)i{xa,Xb) = IVDi - TxxX A\f^s - rx2XB\f^s'^ , 

(\)2{xa,xb) = \yD2 - r22XB^/E~s - r23f{xA,XB)y^s\'^^rid 

(I)3{xb,Xr) = \yD2 - r22XB\f^s - r2ZXR\fE^s'^ • 

The following proposition gives a sufficient condition under 
which Algorithm 1 implements the decoder given in (|4]). 

Algorithm 1 Decoding Algorithm used at D 



for i = 1 to M do 

Xb ^ Si 

Find x\ = arg min {(I)i{xa,xb) + 02(xa, xs)} 

XA e5 

Find x\ = arg min {^i (xa ,xb)} 

Xa g5 

Find Xr = arg min {(I)3{xb, xr)} 

Xr^S 

if (I)i{x\,xb) + (I>2{x\,xb) < (I>i{x\,xb) + (I)3{xb,xr) 

+ \og{SNR) 

then 

m{xB) = + (I)2{x\,xb) 

xa{xb) = x\ 
else 

m{xB) = (I)i{x\,xb) + (t>?>{xB,XR) + \og{SNR) 
xa{xb) = x\ 
end if 
end for 

Find Xb — arg min m{xB) 

XB^S 



15: {xa.Xb) ^ {xa{xb),xb) 



Proposition 1: Algorithm 1 implements the decoder in (|4]), 
if the constants a, 6, c and d are such that the weight matrices 
Wa and Wr are Hurwitz-Radon (H-R) orthogonal, i.e., 
WaWk*+WrWa*=03.[] 

Proof: The decoding metric of the decoder given in (|4]) 
can be written as, 

min {mi {xa.xb) , log (SNR) + m2 (xa, xb)} 

(xa,xb) 

= min {mi (x^, x^) , log (S'A/'i?) + mi (x^, , 

(xa,xb) 

\og{SNR)+m2 {xa,xb)}. 
= min {mi {xa,xb) Aog{SNR) -\- 1713 {xa,xb)}- (9) 

(xa,xb) 

where the metric ms{xA^XB) is given in (|7]), at the top of this 
page. We have, 

m3{xA,XB) = min{||yD^ -HeqxV^lP} 

xrGS 

= min{||yD-Rx/E;f }. (10) 

xrGS 

algorithm exactly similar to Algorithm 1 can be used with the roles of 
A and B interchanged, if Wb and Wr are H-R orthogonal. 



Since Ri is upper triangular, r2i = 0. Also, the entry 
^13 = 0, since Wa and Wr are H-R orthogonal (follows 
from Theorem 2, 1 15 |). Hence, from ([TO]), it follows that. 



(11) 



m2,{xA,XB) = Itoi -riiXA 
+ min \\yD2 -r22XB 

xneS I 

Hence, 

min ms{xA,XB) = min ^i(xa,xb) + min (t)s{xB,XR). 

XA^S XA^S XReS 

From (|9]), the decoding metric can be written as. 



mm < min (0i {xa, xb) + ^2 (^a, xb)) , 

xb^S \xa^S 



min 01 {xa.xb) 

xa^S 



mm ( 



{xB,XR)^\og{SNR) 



In Algorithm 1, inside the for loop, is fixed and the 
operations in lines 3, 4 and 5 entail a complexity order 
0(M). The operations from line 6 to line 12 involve constant 
complexity, independent of M. Hence the complexity order 
for executing the for loop from line 1 to line 13 is (9(M^). 
The operation in line 14 involves a complexity order 0{M). 
Hence the overall decoding complexity order of Algorithm 1 
is (9(M^), which is the same as that of the CFNC scheme 
proposed in p3| . 

Example 2: For the case when a = l, 6=^, c = and 

"1 0" 

the weight matrices are given by, Wa = 



d = 



1 

72' 



Wb 





1 







1 





and Wr 



It can be verified 



that the matrices Wa and Wr are H-R orthogonal, i.e., 
WaWr* + WrWa* = O3. Hence, for this case. Algo- 
rithm 1 can be used to implement the decoder given in (|4]). 

IV. Simulation Results 

Simulation results presented in this section compare the per- 
formance of the proposed PNC scheme with the CFNC scheme 
proposed in p3| . In all the simulation results presented, the 
values of the constants are chosen tobea = l,6 = ^^^ = 
and d = ^ and 4-PSK signal set is used at the nodes. The 
Modulo-4 Latin Square is used at the relay node. 

Fig. |3] shows the SNR Vs. Symbol Error Probability (SEP) 
plots for the case when the variances of all the fading links are 
dB. It can be seen from Fig.[3]that the PNC scheme performs 
better than the CFNC scheme and offers a gain of nearly 3.3 
dB at high SNR. Fig. [4] shows a similar plot for the case when 
the links from A-R and B-R are stronger than the other links, 
i.e., a\j^ = a%j^ = 10 dB, = = ^Id = ^ dB. It 
can be seen from Fig. [4] that for this case, the PNC scheme 
offers a gain of nearly 3 dB at high SNR. Fig. [5] shows the 
plots for the case when the R-D link is stronger than all other 
links, i.e, alj^ = a%j^ = (j% = 



dB. For this case, the PNC scheme offers a gain of about 6.5 
dB at high SNR. Also, it can be verified from the plots that 
the diversity order for the proposed PNC scheme is two. 




SNRindB 



Fig. 3: SNR vs SEP plots for the PNC and CFNC schemes for 
4-PSK signal set for a\j^ = cr|^ = = cr|^ = cr|^ 
dB. 








Fig. 4: SNR vs SEP plots for the PNC and CFNC schemes for 
4-PSK signal set for a\j^ = a%j^ = 10 dB, a\j^ = a 



BD 



^ RD 



dB. 



V. Discussion 



dB, al^ 
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A scheme based on physical layer network coding was 
proposed for the Multiple Access Relay Channel (MARC). 
A novel decoder was proposed for the PNC scheme. The 
conditions which the network coding map / and the constants 
a, 6, c and d should satisfy, for the proposed decoder to offer 
a maximum diversity order two were obtained. It was shown 
that if the constants a, 6, c and d are chosen properly, the 
proposed decoder can be implemented efficiently by a fast 
decoding algorithm. Simulation results presented showed that 
the proposed decoder performs better than the CFNC scheme 
proposed in |13|. The problem of optimizing the choice of 
the constants a, 6, c and d to minimize the error probability 
in addition to ensuring maximum diversity remains open. 



-B-CFNC Scheme: 
-0-PNC Scheme 




Fig. 5: SNR vs SEP plots for the PNC and CFNC schemes 
for 4-PSK signal set for a\j^ = = a\jj = = dB, 



^RD 



10 dB. 



Extending the scheme to a Multiple Access Relay network 
with more than two source nodes and multiple relay nodes is 
a possible direction for future work. 
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APPENDIX - Proof of Theorem 1 

Let H = {Har^ hBR^ ^rd^ ^ad^ ^bd) denote a particular 
realization of the fade coefficients. Throughout the proof, the 
subscript in a probability expression indicates conditioning 
on the fade coefficients. For simplicity of notation, it is 
assumed that the variances of all the fading coefficients are 
one, but the result holds for other values as well. 

Let E denote an error event that the transmitted message 
pair {xa-,xb) is wrongly decoded at D. 

The probability of E conditioned on H given in ([12]), can 
be upper bounded as in ([TSj (eqns. ^V2\ and ([13]) are shown 
at the top of the next page). Ph {/(^A'^b) — 
and Ph {f{x^^x^) ^ f {xa^xb)^ respectively denote 
the probabilities that R transmits the correct and wrong 
network coded symbol during Phase 2, for a given H. 
Also, Ph{{x^,x^) 7^ {xa.xb) 
and Ph {{x^,x^) {xa.xb) 



f{x%x^) = f {xa.xb) 
f{x^,x%) 7^ /{xa^xb) 



respectively denote the probabilities of E given that R 
transmitted the correct and wrong network coded symbol 
for a given H. Ph{E} can be upper bounded as in ([14]), 



where Ph {f{ 



x^^ x^) 



ic'^ } denotes the probability that the 



network coded symbol transmitted by R is x'^ ^ /{xa^xb) 



and Ph{{ 



^ {xa.xb) I f{x 



R 

Ai 



> } is the 



probability of E given that R transmits x'^, for a given H. 
Taking expectation of the terms in (T4\ w.r.t we get ( [TS] ). 

The rest of the proof of Theorem 1 is presented in two 
parts as Lemma 1 and Lemma 2. In Lemma 1, it is shown 



that P{(x2,xg) ^ (xa^xb) 
has a diversity order two. 



f{x^,x§) = /{xa^xb)} 
Lemma 2 shows that 
= x'j^} has a diversity 



P{{x^,X^)y^{xA^XB)\f{x^^ 

order one. Since P{/(f^,x^) = x'^} has a diversity order 
one. Lemma 1 and Lemma 2 together imply that P {E} has 
a diversity order two. 

Lemma 1: The probability P{(f;^,f^) ^ (xajXb) \ 
f{x^^x^) = f (xa^xb)} has a diversity order two. 
Proof: 

Recall that the decoder used at D given in ^ in Section 
II A, involves computation of the metrics mi and m2 defined 
in ([5]) and ([6]). Under the condition that R transmitted the 
correct network coding symbol, a decoding error occurs at 
D only when mi(x^,XB) > mi(x^,x'^) or mi{xA^XB) > 
\og{SNR) +m2(x^,x'^) for some {x'^.xb) 
Hence P{{x^,x^) ^ {xa^xb) \ f{x^,x§) = /{xa^xb)} 
can be upper bounded as in ([16]), which can be upper bounded 
using the union bound as in (17]) (eqns. ([16]) and ( [TT] ) are given 
at the top of the next page). 



Ph{E} = Ph {(^S, ^b) ^ (xa^Xb) I fix^, X^) = f ixA,XB)^ Ph X%) = f ixA,XB)^ 

+ Ph {(^5,^2) ^ (xa^xb) I /(^A' Xb) 7^ / (a^A, a^is)} {f{x^, x^) 7^ / (a^A, ^eb)} (12) 
< Ph {(^5, ^ (xa^xb) I /(icj, ic^) = / (a^A, ics)} + Ph ^{x^^x^) / (a?A, iCis) | fix a, ^b) ^ / (xa^xb)^ Ph {/(^cJ, ^b) ^ / (a^A, a^s)} 



(13) 



Ph {{xa^x^) 7^ (a:A,a;B) | f{x^,x^) = f {xa,xb)] + Yl {(xa^x^) ^ {xa,xb) \ /(^S'^S) = x'r} Ph = ^k} 

^k^/(-A.-i3) 



^A '-^B^"-^/?' 
x'p^^f(xA,XB) 



(14) 



P{£;} < P {(^eJ,:^^) 7^ (xA,a:i3) I /(^2,^b) = / (a^A,a:i3)} + ^ P { (^rj , ^cg) 7^ (xa , a^s) | /(^^^^b) = ^k} ^ {/(^a^^b) = ^k} • (15) 



P |(:e5, :e^) 7^ (a^A, a^s) I /(^J' ^b) = / (xa^xb)^ = P {{mi{xA,XB) > mi{x'^,x'^), {x'^,x^) 7^ (a^A, a^s))} 

U {mi{xA,XB) > log{SNR) + m2(xA' ^s)' (^^A' ^b) / (a^A, a^s)} | /(^A' ^b) = / (a^A, a^is)} (16) 

< ^ P {mi(xA, a^B) > mi(x'^, x'^) I /(^rj, = / (a^A, a^B)} 
/-A'-B)e52 
(=^A'^b)^(^A'^b) 

+ Yl P{miixA,XB)>\ogiSNR)+m2ix'A,x'j,)\f{x^,x^) = f{xA,XB)}. (17) 
(^A '2;'^)es2 

(^A'^b)^(^A'^b) 



7714 (a? A, a^B, a;/?) = I yi^i — ^ AD V^ci a: A — ^bdV^^ a^-B I + |2yi?2 ~ ^-^d a^A — ^bd V^^(i a^B — ^ki? V^a^-R I +log(5A/'P). (18 



' {mi(a;A, a^B) > m2(a?A' ^b) + ^og{SNR) \ /{x^, x^) = f (xa^xb) \ = P < mi{xA,XB) > min 7714 (a^A ' ^^b ' ^k) I /(^A' ^b) = / (a^A, a^B) > 

^ ^ [ ^R^f^^A^^B^ J 



°H \jni{xA,XB) > m^ix^AyXByx'j^) I f{x^,x^) = f{xA,XB)^ 
Ph {\ZD^\^ ^IzD^f >logiSNR) + \zD^ + /lA 



(19) 

< X] P {rni{xA,XB) > m4ix^,x'j^,x'j^) \ f{x^,x^) = f {xa,xb)^ (20) 
^k^-^^^A'^b) 



= Pjf {I + \ZD2\'^ > log{SNR) + \ZD^ + hAD^/E^aAxA + /ibi? V^^A^cb |^ + |2;d2 + hAD^/E^cAxA + hBD\/E~sdAxB + /ikd V^Aa^fil^}. 

(21) 



P/f <^ w < 



log(5Arp) ||x|| 



V2 J • 



(22) 



P {mi{xA,XB) > mi{x'^,x'^) \ /(x^.x^) = /{xa^xb)} 
is equal to the Pair-wise Error Probability (PEP) of a space 
time coded 3x1 collocated MISO system, with the codeword 
difference matrices of the space time code used at the 

aAxA cAxA 
transmitter being of the form bAx b dAx b 

[ f{xA,XB)-f{x'^,x'Q) 

where /S.xa = xa — x'j^^ Ax b = xb — x'^. When 
Axa,Axb 7^ 0, these codeword difference matrices are of 
rank 2, since the restricted codeword difference matrices 
are full rank. When Ax a = and Axb (and also 
Axb = O^Axa ^ 0), the codeword difference matrices 
are full rank, since /{xa^xb) f{xA-,x'^) (otherwise 
the exclusive law given in ^ will be violated). Since the 
codeword difference matrices are full rank, the probability 

P {mi{xA,XB) > mi{x'j^,x'^) \ f{x%x^) = f {xa.xb)} has 

a diversity order two |16) . Let m4 be a metric as 
defined in ([18]), shown at the top of this page. The 



probability P {mi{xA-,XB) > '^2{x'j^^x'^) ^\og{SN R) | 
/(f ^, x^) = / (x^, Xb)} can be written in terms of the 
metrics mi and m4 as in ([19]), which can be upper bounded 
as in ( |2Q| (eqns. ([19]) - ([22]) are shown at the top of this 
page). 



Let Axr 
and Ax b 



f{xA,XB) 

= Xb 



^R-> 



AXA 



XA 



The probability 



Ph {mi{xA,XB) > m4(x'A, x'b, x'r) \ ) = /(xa^xb)} 

can be written in terms of the additive noise zd^ and 



as 



given m 



Let xi = {hADd^XA + hBDbAxB)VE^j ^2 = 
{HadcAxa + hBodAxB + hRDAxR)\/Ws' Also, let zd = 
[zdx ZD2] and X = [xi X2]^ • Then ( [2T] ) can be simplified 



as in ( [22] ), where w = >/2i?e{z^^}, is distributed ac- 
cording to A/'(0, 1). In terms of the Q function, the prob- 
ability Ph[w< "^^^i^jJJ^^ - ^} in ([22]) can be written 



as Q 



logjSNR) W 
L V2||x|| ^ 72 

coefficients. To complete the proof, it suffices to show that 

l logjSNR) ^ 

can be written as, x 

aAxA cAxa 
bAxB dAxB 
AxR 



7211x11 

The vector 



V2 
X 



. Note that ||x|| depends on the fade 
te the proof, it suffices to s 
^ has a diversity order two. 



/Eg [Had ^bd hRo] 



AX 



Since AX AX* is Hermitian, it is unitarily diagonaUz- 
able, i.e, AX AX* = UEU*, where U is unitary and 
"Ai 0" 



I] = 



with Ai > A2. We have ||x|f = 



A2 


SNR hU5]U*h*. Let h = hU = [hi h2 hs]. The vector 
h has the same distribution as that of h, since U is unitary. 

Since the rank of AX is at least one, Ai > 0. We consider 
the two cases where A2 > and A2 = 0. 
Case 1: A2 > 



For this case, upper bounding Q 
which is upper bounded by e 



logjSNR) 
72||x|| 



_V2 

logjSNR) ||x|| 
V2||x|| V2\ 



V2 

, we have 



by 



< e 4 



^i\iSNR\hi\'^+\2SNR\h2\^) 



■ (23) 



Taking expectation w.r.t \hi\ and |/i2|, from (|23]i, we get, 



'logiSNR) 


l|x|| 


. V2||x|| 


" V2 


'\ogiSNR) 


l|x|| 


I V2||x|| ^ 




• A2 =0 





< 



^ has a diversity order two. 



For this case ||x|| = y%^Tm\hi\. Hence 
- \og{SNR) ||x||- 
V2. 
log {SNR) 



Q 



I V2\\k\ 

= Q 



y%jNR\hi 



(24) 



E Q 



V2V\iSNR\hi\ V2 
Let r = \hi\. Taking expectation w.r.t r, from ^24\ , we get, 

log (SNR) 



+ 



f 

P 

J r- 



V2\\x\\ V2. 
log {SNR) 



V^VXiSNRr 



+ 



yXiSNRr 
V2 , 



2re ^ dr 



logjSNR) 
At SNR 



2Q 



log (SNR) VXiSNRr 



V2y^Xi SNRr 



V2 



dr 



f 

J r— 



/ log(SNR) 
XlSNR 



2Q 



log (SNR) ^ VXiSNRr 



V^VXiSNRr 



V2 



In the following, we show that the integrals I i and I 2 
have diversity order two. Note that ^^^^^^^^^ + VXiSNRr ^ 
as a functio n of r, attains the minimum value when r = 
^^XiSNr' minimum value equals ^/2 \og{SNR). 



Since, Q{x) is a decreasing function of we have, 

^ logjSNR) 



^V2VXiSNRr 

we have. 



V2 



h <2Q 



< 



^2\og{SNR)\ [ 

^ Jr=0 

2 / log(SNR) \ 

SNR V J ' 



< Q ^y/2\og{SNR) 

f log(SNR) 



. Hence, 



re ^ dr 



Since for small e ^ can be approximated as 1 — 
X, at high SNR, we have h < sWr ^I^snr • ^^^^^ 



re ^ dr . 



lim 

SNR^oo 

least two. 
Let ro 



■log(^g^^ 



iSNR^ J 



\og{SNR) 



2, Ii has a diversity order at 



XiSNR 

bounded as, /2 < 



iog{SNR) integral I2 can be upper 



Let r' 



logjSNR) 



logjSNR) 
^/WMSNRr 



yXiSNRr 



yXiSNRr 

. As a function of 



rdr. 



V2VX1 SNRr ^ V2 

is monotonically increasing for r > tq. Also, for r > ro, r 

u '^^ ' ^ f / r' ^ J r''' -2 logjSNR) 

can be written m terms or r as, r — ^ 

V^^iSNR \ 

r < 



We have, dr = dr' ^—-^= = , _ , 

V2XiSNR y ^r'2-21og(5Ari?) 

V2X1SNR ' ^2 upper bounded in terms of r' as 



V2XiSNR 

. Since 



Hence /2 < 



Q{rY 
XiSNR 



1 + 



^2 logjSNR) 

, o^ro / , Q{r')r'dr' 

n J ^2 logjSNR) 



y^r'^ -2\og{SNR) 



dr 



1 21 



\iSNrI 



Q(r)r'^ 



y/2iogjSNR) ^r'^ -2\og{SNR) 



dr . 



1 22 

Upper bounding Q{r') hy e ~, I21 can be shown to be 
upper bounded as ^ ^^^2 7 which falls as SNR~^. Upper 

bounding Q{r') by e~^, and using the transformation t = 
r''^ — 2\og{SNR), I22 can be upper bounded as, 

1 r te--2 



/22 < 



XiSNR^ 7o yft^t^2\og{SNR) 
2\og{SNR) 



dt 



POO 

Jo yty 



e 2 



< 



XiSNR^ 7o ^ft^t^2\og{SNR) 
1 r . 1 ' 



e '^dt 



e 2 



XiSNR^ Jo XiSNR^ Jq ^/i 

2^2V2^\og{SNR) 



dt 



XiSNR^ 

where the second inequality above follows from the facts 
that , ^ < 4t and , ^ < 1 for 

^/t+2logjSNR) — Vt ^"^+2 logjSNR) - 

sufficiently large SNR. The last equality follows from 
the fact that ^dt = V2r(l/2) = V^, where 
r(^) is the integral, T{z) = e-H^'^dt. Since, 



log 



2+2v/2^1og(S'Arii:)^ 



XiSNR'^ 



2, 1 22 has a diversity 



lim 

SNR^oo \og{SNR) 
order 2. This completes the proof of Lemma 1. ■ 

Lemma 2: The probability P{(x2,x^) ^ (xa,^;^) | 
f{x%x^)=x'p^^^x'p^ / /{xa^xb)^ has a diversity order 
one. 

Proof: Let 7714 denote the metric as defined in ( [25] ), 
shown at the top of the next page. Under the condition 
that R transmitted the wrong network coded symbol x'^, a 
decoding error occurs at D only when m/^{xA^XB^x'j^ > 
m4,{x'^,x'^,x'^) or m4^{xA,XB,x'^) > mi{x'^,x'^), for 
some {x'^^x'^) / {xa^xb) and x'^ ^ f{x'^^x'^). Hence, 
P{(x2,fg) / (xa,^:^) I =:^'i?} can be upper 

bounded as in ([26]) (eqns. ( [26] ) - ( [29] ) are shown at the top 
of the next page). Using the union bound, from ( [26] ), we get 
(|27j. 



Since 



the 



matrix 



a{xA - x'^) c{xa - x'^) 
h{xB — x'^) d{xB — x'^) 
x' - x" 



has rank at least one for {xa^xb) 7^ 
P {m^{xA,XB, x'^ > m4{x% x'^^.x'^) \ f{x% x%) = x'^} 
has a diversity order at least one. 



Let Axr 

and Axb = xb — 

Ph {m4{xA,XB,x'ji) > mi{xA,XB 



xa 



be written in terms of the additive noise z^i and 



x^ . The probability 

/(^A,^i) = a:^} can 
as 



given in ( [28] ). 

Let xi = {hADdAxA + hBDbAxB)y/E^, X2 = 
{HadcAxa + hBodAxB + hRDAxB)\fWs' Also, let zd = 
[^Di and X = [xi X2Y • Then ( [28] ) can be simplified as 
in ( [29] ), where = >/2Pe{z^^}, is distributed according 
to A/'(0,1). Hence, 



Ph{w< 



log (^TVP) 
V2||x|| 



1 



{||x||2<iog(5Ari?)} 



_ M 
1-Q 



1 



{||x||2>log(5Ari?)}^ 



log(57Vi?) ||x|| 
V2||x|| V2 
log(5A/-i?) ^ ||x|| 
V2 



(30) 



v/2||x|| 

Taking expectation with respect to the fade coefficients in 



30b, P <^ w; < 



P{w< 



{'^ ^ V2||x|| 

log {SNR) 



log(SNR) _ Jjxjl 
V2 



V2\ 



V2 



I can be upper bounded as, 
^<p{\\xf <log{SNR)} 



H:{\\^f>\ogiSNR)} 



log {SNR) 
V2\\x\\ 



The vector x 

/Eg [Had hBD hno] 



can be written 

aAxA cAxa 
bAxB dAxB 
Axr 



M 

V2 



as. 



dH. 



(31) 



AX 



Since AXAX* is Hermitian, it is unitarily diagonalizable, 
i.e, AXAX* = UEU*, where U is unitary and E = 



Ai 0" 
A2 




with Ai > A2. Since the rank of AX is at 



least one, Ai > 0. We have ||x|r = SNR hU5]U*h*. Let 
h = hU = [hi h2 hs]. The vector h has the same distribution 
as that of h, since U is unitary. Hence, we have, 

Plhf < log(SNR)} =p{A,|S,|» + A,|te|» < 



Since, is exponentially distributed, 

P{||xf <log(5iVi?)} < (l-e- 



logjSNR) 
Ai SNR 



)■ 



_ log(SNR) 

At high SNR, 1 — e >^isnr be approximated as 

logjSNR) 
XiSNR ' 



P|||x||^ < log(5'A/'i?)| has a diversity order at 



least one since lim 



\^ XiSNR J 



1 /oArr.N 1. Since (5(x) < 

SNR^oo log (SNR) ^ ^ 

e 2 ^ the integral on the right hand side of ( [31) can be upper 
bounded as, 

log {SNR) ^ ||x||- 
V2. 



L 



Q 



H:{M^>los{SNR)} 

si 

JH:{M\^>\oz{SNR)} 



V2||x|| 

(- ^°l^xP +11-11)' 



dH 



dH. 



(32) 



We consider the following two cases when A2 > and A 2 
0. 

Case 1: A2 > 

For this case, from the integral in ( [32] ), we get, 

log(5iVP) ||x||- 
V2. 



//:{||xf >log(SArK)} 

CX) 



Q 



V2I 



(iff 



< 



SNR 



-\h^\\-\h^\' d\hi\^ d\h2 



(1+A^)(1, ^ 



\2SNR \ 



which falls as at high SNR. 

Cfls^ 2; A2 = For this case, 

\og{SNR) , ||x| 



I Q 

lH:{\\xf>]ogiSNR)} 



V2\ 



V2 



f{H)dH 



< 



\hi\> 



V'log(SWfl) 
,/Ai SNH 



2|^i|e-l'*il'd|^i|. (33) 



The above inequality follows from the fact that for \hi\ > 

y/losjSNR) logjSNR) /WWM 

-vxrmw' ~ vxiSNR\h^\ < -V^('s{SNR). 

Let r = \hi\. From ( [33] ), since (5(^) < the integral / 



'm4ixA,XB,XR) = lym - hAD^/E^a XA - hnoVE^b xb\ -\- \yD2 - J^adV^^c xa - hBD^/E^d xb - hnoVE^XRl -\-log{SNR). (25) 



P {(xa^^b) 7^ ixA,XB) I fix^,x^) = x'^} = ^ P {{m4ixA,XB,x'^) > m^ix'J^ , x^ , x^) , x'^ ^ f{x'^,x'^)} 

U {m4ixA,XB,x'p>) > mi(a?'^, a^'j^)} | /{x^, x^) = x'^^^ (26) 

< ^ P {m4ixA,XB,x'^) > m4ix'X,x^,x^) \ fix^,x^) = x'^^ 

(x'J,,x-)es^ x'^es, 

{x'^,x'^)^{xa,xb) a:'^^/(^A'^B) 

+ P{m4(xA, iCB, ic'^) > mi(a?'^, a?'^) | /(^rj, f^) = } • (^7) 



P/f |log(5Ari?) + m4{xA,XB,x'j^) > mi(x'^, x'j^) | /(x^, ic^) = x'^^ j 



Ph{\zd^\ + l^^i^sl +log(5ArP) > \zD^ -^hADVEsaAxA + hBDVEsbAxBl + -^hADVEscAxA -^hBDVEsdAxB -^hRoAxRl^}. 

(28) 



P /oj? / * ^ \ . logOSTVP) \ r ^ log jSNR) ||x|| 

Pi^ < 2Re < zo > < ^ X > = P/f < w < — — 

I I l|x|| J - ||x|| "J \ - y2||x|| V2 



(29) 



can be upper bounded as, 

XiSNR\ 



I < 



^yXlSNR\og(SNR) 
A + XiSNR 



y{SNR) 



log(SNR) log(SNR)XiSNR 
g 4 + 16 + 4XiSNR 
^ V ^ 

(34) 



Let ro, /ci, /c2 and /cs be defined as shown in ( |34| ). From 
([34]), the upper bound on / can be written as. 



r>ro «^ ^^^0 



/ < Y / 2(r - fc2)e-'=i('-'==)'dr + k2 kg [ 

^ J r>ro «^ 

Since X,>^^ e-'=i('^-'=^)'dr = [(ro - k2)V2h] 

I < -^e-'=^('-°-'==)' + k2 ks.f^Q [(ro - fc2)v^ 
zki V ^1 

Since Q [(ro - k2)V2h] < 1, e-^^^^o-^^)' < ^nd ks 
can be approximated as one at high SNR^ substituting for 
ro, ki,k2 and ks, we get, 

1 , v'7rAi5Ari?log(5Ari?) 



2 (1 + ^) (4 ^ ^^sNR) Jl + Ai^ 



Since the above upper bound on / falls as SNR~^ at high 
SNR, / has a diversity order at least 1. This completes the 
proof of Lemma 2. ■ 



